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摘要 采用 有 限 元 法 和 热 介 


怨 张 量 模拟 研究 了 非 


波 传播 和 功率 沉积 的 影响 。 模 拟 结果 显示 : 


次 吸收 假设 , 以 及 一 次 吸收 假设 对 EAST 四 电流 带 天 线 少数 离子 加 热 


在 其 他 参数 相同 的 条 件 下 ， 当 天 线 电 流 相位 为 (0 0 x m) AO x x 0) 时 ， 两 种 模型 计 
算 的 天 线 附 近 电 场 分 布 基本 一 致 ， 当 天 线 电流 相位 为 (0 x 0 z) 时 ， 两 个 模型 的 电场 模拟 结果 出 现 显著 差别 。 环 向 传播 多 圈 的 


离子 回旋 波 是 导致 结果 差异 的 主要 原因 。 环 形 聚 变 装置 中 ， 当 波 的 平行 波 数 较 大 时 ， 一 次 吸收 假设 的 模拟 结果 偏 小 。 对 于 环 


形 聚 变 装置 ， 一 次 吸收 假设 的 适用 
关键 词 ”离子 回旋 波 ， 少 数 五 氢 离 子 加 热 ， 功 率 沉 积 ， 热 介 电 张 量 ， 完 美 匹配 层 


中 图 分 类 号 ”0532 


性 需要 进一步 系统 的 研究 。 


DOI: 


离子 回旋 共振 加 热 (CRH) 是 大 多 数 仿 星 器 和 托 卡 马 死 装置 中 使 用 的 主要 辅助 加 热 方法 之 一 03。 对 于 较 


大 或 开放 的 等 离子 体 域 ，ICRH 程序 通常 会 使 月 


日 吸收 边界 条 件 截 断 求解 域 。 这 种 截断 方法 也 被 称 为 一 次 吸 


收 假设 。 常见 的 吸收 边界 条 件 有 : 辐射 边界 由 、 完 美 匹 配 层 ML) 向 和 波动 方程 的 单 向 近似 方法 加 。Lamalle 


等 人 四 开发 了 初步 的 ICRH 耦合 程序 , 程序 中 使 
响 ，Jacquot 等 人 中 开发 了 ICRH [IPEA EU, RUFI T PML 代 蔡 芯 部 等 离子 体 对 离子 


pu 


海水 介质 代替 了 等 离子 体 并 利用 PML 减少 了 反射 波 的 影 


| 


旋 波 的 


吸收 ; Shiraiwa 等 人 四 开发 了 WEST 装置 的 ICRH 加 热 程序 , 程序 还 原 了 WEST 装置 的 天 线 和 部 分 真空 室 
结构 ， 并 在 真空 室 的 两 端 使 用 了 辐射 边界 。 


在 聚变 装置 中 ， 一 次 吸收 假设 意味 着 波 在 绕 真 空 室 传 播 一 圈 之 前 将 会 被 等 离子 体 完全 吸收 ,但 该 假设 


在 聚变 装置 中 的 适用 性 是 缺少 系统 研究 。 本 课题 基于 COMSOL 软件 ， 开 发 了 非 一 次 吸收 假设 的 EAST 全 
真空 室 模 型 和 一 次 吸收 假设 的 EAST 1/4 真空 室 模型 ， 并 在 (0r0m，(00rm，(0rr0) 三 种 天 线 电流 相位 
下 对 比 研究 了 两 种 模型 在 D(D) 等 离子 体 少数 离子 加 热 时 离子 回旋 波 的 传播 和 功率 沉积 。 


1 物理 模型 和 方程 


11 EAST 三 维 全 真空 室 ICRH 程序 


EAST 三 维 全 真空 室 域 离子 回旋 加 热 程 序 的 结构 示意 图 和 EAST 真空 室 二 维 截面 图 如 图 1 所 示 。 其 
空 室 结构 参照 了 目前 EAST 正在 使 用 的 带 有 下 钨 直角 偏 滤器 的 第 一 壁 结构 。 
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离子 回旋 天 线 放置 在 x 


AR 


正方 向 约 2.4 m 位 置 处 。 模 型 中 未 使 用 一 次 吸收 假设 。 


图 1 EAST 三 维 全 真空 宝 域 离子 回旋 加 热 程序 的 结构 示意 图 


Fig. 1. Schematic diagram of the structure of the EAST three-dimensional full-vacuum chamber ion cyclotron heating program 
离子 回旋 波 天 线 结构 参照 了 EASTI 窗口 四 电流 带 ICRH 天 线 ， 如 图 2 所 示 。 图 2(a) 中 ， 分 别 用 红色 、 
灰色 和 橙色 标记 了 电流 带 、 法 拉 第 屏 和 天 线 盒 。 极 向 曲率 半径 约 760 mm， 环 向 曲率 半径 约 1980 mm. X 
线 盒 中 间 安 装 了 三 块 隔 板 以 避免 电流 带 与 电流 带 间 的 耦合 。 天 线 盒 的 背面 有 四 个 传输 线 (VTL) 的 接口 ， 四 
根 传 输 线 分 别 连接 四 个 电流 带 (图 2b))。 电 流 带 的 折 角 处 统一 使 用 半径 为 30 mm 的 圆 角 设计 ， 为 了 增强 电 
流 带 散 热 ， 电 流 带 采用 了 炙 空 设计 。 天 线 和 天 线 箱 的 材料 设计 为 316L 钢 ， 设 置 为 阻抗 边界 ， 有 热 损耗 。 
目前 天 线 结构 中 暂 未 包含 限制 器 结构 。 


(b) 


septum FS antenna box return conductor 


图 2 带 有 FS 的 天 线 视图 (a) 与 四 电流 带 、 返 回 导体 及 VTL 的 视图 (b) 
Fig. 2. View of the antenna with FS (a) and View of the four current straps, return conductor and VTLs(b) 
1.2 EAST 三 维 1/4 AZ Z ICRH 程序 
图 3 为 EAST 三 维 1/4 ATZ ICRH 程序 的 结构 示意 图 ， 其 中 灰色 区 域 代表 真空 室 ， 其 环 向 角 覆 盖 范 
围 为 -r/4 至 mw/4。 真 空 室 的 两 端 为 PML( 绿 色 区 域 )。PML 可 以 无 反射 地 吸收 进入 该 区 域 的 电磁 波 以 模拟 环 
向 无 限 长 的 真空 室 结构 ， 即 一 次 吸收 。 天 线 结构 与 EAST 全 真空 室 ICRH 程序 中 相同 。 


图 3 EAST 三 维 全 真空 室 域 离子 回旋 加 热 程 序 的 结构 示意 图 
Fig. 3. Schematic diagram of the structure of the EAST three-dimensional 1/4 vacuum chamber ICRH program 


13 ”等 离子 体 参 数 设置 

EAST 四 电流 带 天 线 与 等 离子 体 相 互 作 用 时 ， 天 线 私 有 区 (图 1 蓝 色 部 分 ) 为 真空 区 域 ， 真 空 室内 充满 
等 离子 体 ( 图 1 灰色 部 分 )。 等 离子 体 密度 和 温度 采用 EAST#107399 炮 实验 数据 (图 和， 磁场 使 用 双 零 位 型 ， 
采用 低 杂 波 (CHW) 加 热 、 电子 回旋 共振 加 热 CCRH) 和 离子 回旋 共振 加 热 (CRH) 等 手段 加 热 等 离子 体 (图 4)。 
采用 EFIT 和 SOLPS 程序 重建 #107399 炮 8s 时 的 磁场 分 布 (如 图 5(a))、 密 度 分 布 (如 图 5(b))、 电 子 温度 分 布 
(如 图 5(c)) 和 离子 温度 分 布 (如 图 5(d))。 在 氧气 等 离子 体 中 采用 少数 氧 离子 基 频 加 热 ， 氧 浓度 为 10%。 天 线 
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4 EAST 脉冲 107399 期 间 参数 随 时 间 的 变化 
Fig. 4. Temporal behavior of parameters during the EAST pulse 107399 
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5 ” (a) 磁场 重建 结果 ; (b) 
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包子 密度 重建 结果 ，(c) 电 子 温度 重建 结果 ; 


(d) 离 子 温度 重建 结果 


Fig5. (a) Magnetic field reconstruction results; (b) Electron density reconstruction results; (c) Electron temperature reconstruction 


results; (d) Ion temperature reconstruction results 


14 等 离子 体 设 置 
磁化 等 离子 体 中 的 波动 方程 为 : 
v(v 


式 中 : 
子 体 的 等 效 相对 介 电 张 量 表示 为 外 
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式 中 :n 为 谐 波 级 数 ，a ”为 不 同 的 粒子 , I Q AIO, 为 n 阶 修正 贝 塞 尔 函 数 及 其 导数 ; Z(6 ) 
RI Z'(£) 为 等 离子 体 色散 函数 及 其 导数 ，6 ,= o nO, /ves 为 色散 函数 的 阶 参数 。 
EAST 四 电流 带 天 线 主要 使 用 (0 0 m m). (0 x x 0) 和 (0 x 0 7) 三 种 电流 相位 加 热 等 离子 体 。 它们 的 功率 谱 


峰值 分 别 为 k | = 5rad/m.k | =8rad/m 和 k | = 13 rad/m。 当 等 离子 体 密 度 较 高 时 ,只 有 快 波 (FW) 
竺 等 离子 体 中 传播 ， 快 波 重 直人 磁场 的 波 矢 分 量 kK 的 大 小 可 以 根据 色散 公式 计算 得 到 钙 ; 
o ko KK, 
i S an = ke 十 KK. F (9) 
4;=0 


热 介 电 张 量 考虑 四 阶 有 限 拉 莫 尔 半 径 近 似 , 即 公 式 (3)-(8) 中 n=[-4, 4]。 采 用 文献 5 中 介绍 的 张 量 旋转 
方法 ， 热 介 电 张 量 中 的 磁 季 方向 可 以 从 z 轴 旋转 到 磁力 线 方向 。 结合 三 维 磁 场 分 布 ， 得 到 托 卡 马克 真空 室 
中 的 三 维 热 等 离子 体 介 电 张 量 2 to 波 在 等 离子 体 中 的 功率 沉积 大 小 为 : 


Py, = | $06 ji (10) 


2 结果 与 分 析 


2.1 电场 对 比分 析 
本 节 中 ， 使 用 case 1 代表 EAST 三 维 1/4 真空 室 ICRH 程序 ，case 2 代表 EAST 全 真空 室 ICRH 程序 。 

图 6- 图 11 对 比分 析 了 三 种 相位 条 件 下 两 种 程序 模拟 的 离子 回旋 波 在 等 离子 体 中 的 传播 。 图 6(a)- 图 
11(a) 为 离子 回旋 波 x-y 截面 电场 分 布 图 ， 图 6(b)- 图 11(b) 为 离子 回旋 波 x-z 截面 电场 分 布 图 。 模 拟 结果 表 
明 ， 电 流 相位 为 (0 0 zx) 和 (0 x 0 DN, WA cases 模拟 的 天 线 附近 电场 分 布 基本 一 致 ， 电 流 相 位 为 (0 x 0 m) 
时 , 两 个 cases 模拟 的 天 线 附 近 电 场 分 布 差异 增 大 。 可 能 的 原因 是 平行 波 数 越 大 , 波 传播 方向 与 磁场 方向 ( 环 
向 ) 的 夹 角 越 小 。 夹 角 越 小 ， 离 子 回 旋 波 沿 环 向 传播 距离 越 长 ， 有 些 离子 回旋 波 不 满足 一 次 吸收 假设 可 以 围 
绕 真 空 室 传播 多 圈 ( 图 11(a))， 天 线 附近 电场 是 发 射 波 和 部 分 绕 环 向 转 一 圈 或 多 圈 波 的 县 加 ， 因 此 case 1 和 
case 2 的 模拟 电场 分 布 不 同 。 
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图 6 Case 1 中 EAST 天 线 电流 相位 为 (0 0 x zn) 时 的 左旋 波 电 场 分 布 (a) x-y 截面 ; (b) x-z 截面 


Fig. 6. The distribution of the right-handed electric field of the EAST antenna in case 1 when the current phase is (0 0 2 x) (a) x-y 
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图 7 Case 2 中 EAST 天 线 电流 相位 为 (0 0 c o RED Ze BS HIP (a) x-y 截面 ; (b) x-z 截面 
Fig. 7. The distribution of the right-handed electric field of the EAST antenna in case 2 when the current phase is (0 0 x x) (a) x-y 
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8 Case 1 中 EAST 天 线 电 流 相位 为 (0rr0) 时 的 左旋 波 电场 分 布 (a) x-y 截面 ; (b) x-z BETH 
Fig. 8. The distribution of the right-handed electric field of the EAST antenna in case 1 when the current phase is (0 x x 0) (a) x-y 
sections; (b) x-z sections 
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9 Case 2 F EAST 天 线 电流 相位 为 (0 x x 0) 时 的 左旋 波 电 场 分 布 (a) x-y 截面 ; (b) x-z ATE 


Fig. 9. The distribution of the right-handed electric field of the EAST antenna in case 2 when the current phase is (0 x x 0) (a) x-y 
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图 10 Case 1 中 EAST 天 线 电 流 相 位 为 (0 r0m) 时 的 左旋 波 电场 分 布 (a) x-y 截面 ; (b) x-z 截面 


Fig. 10. The distribution of the right-handed electric field of the EAST antenna in case 1 when the current phase is (0 2 0 x) (a) x-y 
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图 11 Case 2 中 EAST 天 线 电 流 相位 为 (0 x 0 7 时 的 左旋 波 电场 分 布 (a) xy 截面 ，(b) x-z ARTE 


Fig. 11. The distribution of the right-handed electric field of the EAST antenna in case 2 when the current phase is (0 ™ 0 1) (a) 


x-y sections; (b) x-z sections 
2.22 ”功率 沉积 分 析 
表 1 展示 了 在 case 1 和 case 2 等 离子 体 总 功率 沉积 的 计算 结果 。 由 表 1 可 知 ，case 1 和 case 2 在 相同 
电流 相位 下 功率 沉积 相对 误差 不 超过 2%。 该 现象 说 明 ， 在 平行 波 数 k 1 较 大 时 一 次 吸收 假设 会 使 电场 的 
预测 变 差 ， 但 对 总 功率 沉积 的 预测 基本 符合 。 此 外 ， 对 比 不 同 电流 相位 的 总 功率 沉积 发 现 ， 天 线 与 等 离子 
体 的 耦合 功率 随 K 1 的 增 大 而 减 小 。 电 流 相位 对 耦合 功率 的 影响 与 EAST 上 的 实验 结果 趋势 上 是 一 致 的 [1]。 
表 1 Case 1 和 case 2 不 同 电流 相位 下 的 总 功率 沉积 


Table 1. Total power deposition at different current phases in Case 1 and Case 2 


Case (0027) (0 x 20) (0 20 x) 
Case | 0.81 MW 0.75 MW 0.65 MW 
Case 2 0.82 MW 0.74 MW 0.64 MW 


表 2 展示 了 在 不 同 天 线 电流 相位 下 ， 真 空 室内 电子 的 功率 沉积 份额 。 由 表 1 AYA, case 1 和 case 2 的 
电子 功率 沉积 都 随 着 增加 而 增 大 。 可 能 的 原因 是 随 着 k “| 的 增 大 ， 电 子 对 离子 回旋 波 的 阻尼 效应 增强 ， 并 
且 更 大 的 k “| 可 能 意味 着 离子 回旋 波 的 功率 沉积 区 域 更 大 。 相 同 天 线 电 流 相位 下 ，case 1 和 case 2 中 电子 
功率 沉积 份额 的 差别 也 随 着 增加 而 增 大 ， 这 说 明 一 次 吸收 假设 会 影响 粒子 的 功率 沉积 份额 。 

表 2 Case 1 和 case 2 不 同 电流 相位 下 的 电子 功率 沉积 份额 


Table 2. The proportion of electron power deposition at different current phases in case 1 and 2 


Case (0027) (0 x 20) (0 20 x) 
Case | 0% 2% 9.4% 
Case 2 0% 3.5% 13.8% 


ERREIAN x 07) 时 , 图 12 和 图 13 对 比 case 1 和 case 2 的 电子 和 氢 离 子 在 x-y 截面 的 功率 沉 
积分 布 。 所 离子 的 功率 沉积 主要 分 布 在 所 离子 回旋 共振 层 附近 ， 电 子 的 功率 沉积 则 更 集中 于 等 离子 体 的 高 
密 高 温 区 域 。 对 比 相 同 的 1/4 真空 室 区 域 ，case 2 中 电子 和 氧 离 子 的 功率 沉积 密度 明显 较 高 。 可 能 的 原因 
是 case 2 中 不 满足 一 次 吸收 假设 的 离子 回旋 波 产生 了 额外 的 功率 沉积 。 结 果 表 明 ， 环 形 装置 会 使 部 分 天 线 
发 射 的 离子 回旋 波 回 到 天 线 附 近 ， 在 部 分 条 件 下 一 次 吸收 假设 的 模拟 结果 会 出 现 较 大 误差 。 
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图 12 Case 1 中 天 线 电 流 相 位 为 (0 x 0 7) 时 的 功率 沉积 分 布 (a) 电子 ; (5) A 


Fig. 12. The power deposition distribution in case 1 when the antenna current phase is (0 x 0 7). (a) Electron; (b) H ion. 
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图 13 Case 2 中 天 线 电 流 相 位 为 (0 m 0 mI IAS DUERA) a) 电子 ; (b) 氧 离子 


Fig. 13. The power deposition distribution in case 2 when the antenna current phase is (0 n 0 x). (a) Electron; (b) H ion 


3 结论 


采用 EAST 实验 参数 ， 本 课题 开发 了 一 次 吸收 假设 的 EAST = 46 1/4 真空 室 ICRH 程序 和 非 一 次 吸收 
假设 的 三 维 全 真空 室 ICRH 程序 ,使 用 两 种 程序 模拟 了 D(H) 等 离子 体 少数 离子 加 热 时 离子 回旋 波 的 传播 和 
功率 沉积 。 

模拟 结果 显示 : 非 一 次 吸收 假设 中 包含 了 离子 回旋 波 环 向 传播 一 圈 或 多 圈 后 对 电场 分 布 的 影响 。 环 向 
传播 一 圈 或 多 圈 的 离子 回旋 波 占 比 随 平行 波 数 的 增 大 而 增 大 。 不 满足 一 次 吸收 假设 的 离子 回旋 波 会 对 电场 
分 布 和 粒子 功率 沉积 分 布 产 生 明 显影 响 ， 但 对 总 功率 沉积 的 影响 较 小 。 
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Abstract Background: The propagation and power deposition process of ion cyclotron waves in plasma is the key to ICRH. Purpose: 
The purpose of this study is to evaluate the effect of single absorption assumption and non-single absorption assumption on ion cyclotron 
wave coupling and heating in the EAST. Methods: Under different antenna current phases, the simulation results of the EAST 3D full 
vacuum chamber ICRH program and the 1/4 vacuum chamber ICRH program were compared. Results: When the antenna current phase is 
(0 x 0 x), there is a significant difference between the simulation results of the two models. Conclusion: The ion cyclotron wave 
propagating multiple times in the toroidal direction is the main reason for the difference in results. For fusion devices, the applicability of the 
single absorption assumption needs further systematic research. 
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